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ABSTRACT 

We present radio and X-ray observations, as well as optical light curves, for a subset of 26 BL Lac 
candidates from the Sloan Digital Sky Survey (SDSS) lacking strong radio emission and with z < 2.2. 
Half of these 26 objects are shown to be stars, galaxies, or absorbed quasars. We conclude that the 
other 13 objects are Active Galactic Nuclei (AGN) with abnormally weak emission features; ten of 
those 13 arc definitively radio-quiet, and, for those with available optical light curves, their level of 
optical flux variability is consistent with radio-quiet quasars. We cannot exclude the possibility that 
some of these 13 AGN lie on the extremely radio-faint tail of the BL Lac distribution, but our study 
generally supports the notion that all BL Lac objects are radio-loud. These radio-quiet AGN appear 
to have intrinsically weak or absent broad emission line regions, and, based on their X-ray properties, 
we argue that some are low-redshift analogs to weak line quasars (WLQs). SDSS BL Lac searches are 
so far the only systematic surveys of the SDSS database capable of recovering such exotic low-redshift 
WLQs. There are 71 more z < 2.2 radio-quiet BL Lac candidates already identified in the SDSS 
not considered here, and many of those might be best unified with WLQs as well. Future studies 
combining low- and high-redshift WLQ samples will yield new insight on our understanding of the 
structure and formation of AGN broad emission line regions. 
Subject headings: BL Lacertae objects: general — quasars: emission lines 



1. INTRODUCTION 

Strong optical and ultraviolet (UV) emission lines are 
among the most prominent and often defining charac- 
teristics of Active Galactic Nuclei (AGN). The rare in- 
stances of AGN lacking strong emission features are 
typically classified as BL Lacertae objects. In addi- 
tion to their featureless optical spectra, BL Lac objects 
are strong radio, X-ray and gamma-ray emitters, they 
are highly polarized, and they display strong and rapid 
flux variability across the entire electromagnetic spec - 
trum (e.g., see iKollgaardl 119941 : iPerlman et~aH 120011 ). 
In the standard unification paradigm, BL Lac objects 
are explained as low-luminosity radio galaxies with a 
relativistic jet pointed toward the observer ( e.g., see 
IBlandford fc Reesl[l978T: lUrrv fc Padovani|[l995l ). In this 
scenario, BL Lac objects' emission lines are diluted by 
the Doppler boosted relativistic jet, and the BL Lac phe- 
nomenon should go hand in hand with strong radio emis- 
sion. Observations have indeed established that radio- 
quiet BL Lac objects must be extraordinarily rare if they 
exist at all (e.g., se e Uannuzi et alj|1993t ILondish et al.l 
I2004D . For example, IStocke et al.l (| 1990) found no exam- 
ples of radio-quiet BL Lac objects in the X-ray selected 
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BL Lac sample from the Einstein O bservatory Extende d 
Medium- Sensitivity Survey (EMSS. IStocke et al.|[l99ll ). 

It is possible, however, that venerable BL Lac samples 
such as the EMSS were merely too small (~40 objects) to 
reveal especially rare subsets of radio-shy BL Lac objects. 
With large-scale multiwavelength surveys and inclusive 
selection approaches, the number of known BL Lac ob- 
jects has grown over th e past 10-15 years, from a c ou- 
ple hun dred cataloged i n Padova ni fc GiommJI (|1995bl) to 
~10 3 in lMassaro et al.l (|2009D . Quasar samples have also 
become extraordinarily large; for example, the Sloan Dig- 
ital Sky Survey (SDSS, lYork et al.ll2000D Data Release 
Seven (DR7) Quasar Catalog con tains ~10 5 quasars wit h 
reliable spectroscopic redshifts (Schneid er et al.| [2010). 
These quasar samples are now large enough to reveal 
interesting examples of rare quasar populations, such 
as radio-quiet objects with weak or absent spectral fea- 
tures, in num bers comparabl e to ea rly BL Lac samples. 
For example. iCollinge et ail (|2005l hereafter C05) and 
lAnderson et al.l (|2007l . hereafter A07) discovered about 
two dozen BL Lac c andid ates lacking strong radio emis- 
sion. iPlotkin e t al. (2010J, hereafter P10) present a sam- 
ple of 723 optically selected BL Lac objects uniformly 
selected from the SDSS that includes 86 objects with ra- 
dio fluxes or limits firmly placing them in the radio-quiet 
regime. 

This population of weak-featured radio-quiet objects 
recovered by SDSS BL Lac searches is intriguing, and it is 
puzzling how nature can create weak-featured AGN that 
are not strong radio emitters. If these objects' emission 
lines are not simply diluted by beamed emission, then 
some other mechanism must be invoked to explain their 
apparently anemic broad emission line regions (BELRs). 

Around one-third of the weak-featured radio-quiet 
AGN recovered by SDSS BL Lac searches have z > 
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2.2. These high-redshift objects pass the phenomeno- 
logical criteria to alternatively be classified as weak line 
quasars (WLQs) - a rare class of high-redshift AGN dis- 
covered by the S PSS with weak or absent UV emis- 
sion features (seelFan et a.1.1 1 1 999t I Anderson et al.l 2001; 
IShemmer et~aLll2006ll2009HDiamond-Stanic et alj |2009) . 
There is significant overlap in the optical properties 
of WLQs and BL Lac objects, but WLQs tend to be 
weaker radio and X-ray emitters and they show less flux 
variability and polarization. WLQs appear to compose 
an exotic population of fe atureless AGN distinc t from 
BL Lac objects (e.g., see lShemmer_eiaL| HM 120091: 
iDiamond-Stanic et al.ll2009h . That is, they seem to have 
intrinsically weak BELRs, and their emission lines are 
unlikely to be simply diluted by a beamed relativistic 
jet. 

WLQ redshifts can often be determined only from the 
onset of the Lya forest in their SDSS spectra. For 
this reason, studies that explicitly search for WLQs in 
the SDSS database only target high-redshift sources, 
and ~80 high-redshift WLQs have been discovered by 
the SDSS to date. SDSS BL Lac searches are sensi- 
tive not only to high-redshift WLQs but also to lower- 
rcdshift analogs, if such objects existQ Perhaps some of 
the lower- redshift weak-featured radio-quiet AGN from 
SDSS BL Lac searches are actually lower-redshift analogs 
to high-redshift WLQs. 

In this paper we discuss the multiwavelength proper- 
ties of the population of 26 lower-redshift ( z < 2 .2) radio- 
quiet BL Lac candidates recovered by t he IC05I and IA071 
surveys. We briefly discuss objects from lPloT as well, but 
we defer detailed description of their properties for an- 
other paper. Throughout, we refer to these objects as 
"radio-quiet BL Lac candidates" to reflect the method 
by which they were discovered, but we will show that 
not all objects should be unified with normal radio-loud 
BL Lac objects. In Section [5] we present the sample of 
26 low-redshift radio-quiet BL Lac candidates. We then 
present new radio and X-ray observations of a selected 
subset in Sections 13.11 and 13.21 respectively. In Section 
13.31 we discuss their optical variability, and we compare 
them to radio-loud BL Lac objects as well as to normal 
radio-quiet quasars. Contaminants are identified in Sec- 
tionUJ in Section[5]we discuss potential physical interpre- 
tations for their weak-featured spectra, we compare them 
to other lineless AGN in the literature. Our main con- 
clusions are summarized in Section |6] We adopt the fol- 
lowing cosmology: Hg = 71 km s _1 Mpc -1 , Q, m = 0.27, 
and A = 0.73. 

2. THE PARENT RADIO-QUIET BL LAC SAMPLE 

In Table [T] we present all 26 low-redshift ( z < 2.2) 
radio - quiet BL Lac candidates recovered by IC05I and 
IA07L Throughout the text we identify objects that 
should be removed from this list (also see Section 2]), 
but we include all 26 objects here for completeness. 
Sources with z > 2.2 are not included in Table [1] 
because WLQs have been discussed elsewhere (e.g., 
IShemmer et all 120091 : Ibiamond-Stanic et al.ll2009ft . 

We refer the reader to IC05I and IA07I for detailed de- 



7 BL Lac surveys do not recover all known WLQs because BL Lac 
searches typically impose more stringent criteria on spectral feature 
strengths. 



scriptions of their BL Lac selection. Briefly, both sur- 
veys require objects with optical SDSS spectra to not 
show any emission features with rest-frame equivalent 
widths (REW) stronger than 5 A, and all spectra must 
also show a Ca II H/K d epression smaller than 40% (see 
lLandt et al.ll2002f) . IA07I require the additional constraint 
that all BL Lac candidates must match within 60" to 
an X-ray source in the ROSAT All Sky Survey (RASS, 
iVoges et al.1 119991 l2000fl , with no constraints on their 
radio emission. IC05I do not explicitly require radio or 
X-ray emission for inclusion; post-selection, their sam- 
ple is correlated to RASS in the X-ray and to the Faint 
Images of the Rad io Sky at Twenty-cm survey (FIRST, 
Becker et al.lll995h and to the NRAO VLA Sky Survey 
(NVSS. ICondon et al.lll998fl in the radio. Both radio sur- 
veys were performed with the Very Large Arraj0 (VLA) 
at 1 .4 GHz. The 20 radio-quiet BL Lac candidates from 
lG05l are not detected in the radio by FIRST/NVSS (or 
are not in the FIRST footprint), or they have ra dio d etec- 
tions but are radio-quiet. Seven objects fro m IA07I lack 
radio detections in FIRST/NVSS, and one IA07I object 
was also identified as a radio-quiet BL Lac candidate by 

IcofiL 

The parameters lis ted in Table \H are tak en from the 
SDSS DR7 database (jAbazaiian et al.ll2009D . We adopt 
the spectroscopic redshift from the literature, unless we 
can identify a better redshift from visual examination of 
their DR7 spectral reductions. Of the 26 objects, 13 have 
spectroscopic redshifts identified from weak spectral fea- 
tures. Five of the 13 redshifts are tentative (i.e., they are 
less certain because their weak emission line identifica- 
tions are less secure or they show only one emission line 
that we assume to be Mg II), which we treat as exact. 
One of these tentative redshifts is actually a lower limit 
derived from possible intervening Mg II absorption. All 
objects in Tabled] have redshift upper limits z < 2.2, or 
else the Lyman forest would be observed in their SDSS 
spectra. We also list morphology and proper motion in- 
formation in Table [T] to help assess if any objects lacking 
redshifts might be stars. Proper motions are taken from 
the S DSS+USNO-B proper motion catalog (jMunn et al.l 
I2004D . as listed in the SDSS DR7 database. 

Table[T] includes broad-band radio to optical (a ro ) and 
optical to X-ray {a ox ) spectral indices for each entry0 
The a ro values in Table Q] are based on FIRST detec- 
tions (or limits) and SDSS, and the a ox values are de- 
rived from RASS and SDSS. We assume objects with 
unknown redshifts have z = 0.3 and redshift lower limits 
are exact, but we note that a ro and a ox are relatively in- 
sensitive to the exact redshift used. Blank entries fo r a ro 
indi cate the source is outside the FIRST footprint. IC05I 
and IA07I use slightly different parameters for estimating 
broad band spectral indices, so we recalculate them here. 
If contamination from the host galaxy to the SDSS opti- 
cal spectrum is observed, then we decompose the SDSS 

8 The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 

9 The broad-band spectral index, for V2 > v\, is de- 
fined as a v , V2 = —\ofr{L V2 /L vl )/\og{v2/v\). Here, a ro = 
-log(L o /L r )/5.08 and a ox = -log(L x /L o )/2.60, where L r , L , 
and L x are the specific luminosities (per unit frequency) at rest- 
frames 5 GHz, 5000 A, and 1 keV respectively jTananbau m et al.l 

115751 ; IStocke et al.|[T985T) . 



Radio-quiet Quasars with Weak Emission Lines 



TABLE 1 
Radio-Quiet BL Lao Candidates 



SDSS Name 


RA 


Dec 


z sp<!c a 


M b 


r c 


logKL ) d 


<*£••" 


a ox SS 


Class. 


Comments 


(.72000) 


(J2000) 


(J2000) 






(mag) 


(erg s" 1 ) 










012155.87-102037.1 


20.48281 


-10.34366 


0.470 


7.6 


19.52 


44.44 


<0.234 


>1.223 


U 


C05.VLA 


013408.95+003102.4 


23.53731 


0.51735 




147.4 


19.90 




<0.266 


>1.152 


S 


C05,VLA,LC 


020137.70+002535.0 


30.40709 


0.42641 




175.9 


19.50 




<0.235 


>1.238 


s 


C05.VLA 


024156.37+004351.6 


40.48490 


0.73100 


0.990 


3.6 


19.59 


45.13 


<0.217 


>1.074 


u 


C05,VLA,LC 


024157.36+000944.1 


40.48902 


0.16226 


0.790? 


54.1 


20.54 


44.64 


0.249 


0.863 


RLBL 


C05,A07,VLA,LC 


025046.47-005449.0 


42.69366 


-0.91361 




10.9 


19.99 




<0.283 


>0.931 


S 


C05,VLA,LC 


025612.47-001057.8 


44.05197 


-0.18273 






20.36 




<0.288 


>0.876 


u? 


C05,VLA,LC 


075551.44+352549.8 


118.96436 


35.43050 




9.1 


19.06 




<0.208 


1.110 


X 


A07 


090133.42+031412.5 


135.38927 


3.23681 


0.459 


1.6 


18.96 


44.56 


<0.203 


>1.224 


u 


C05 


100847.01+114946.9 


152.19591 


11.82972 


0.261 


11.3 


19.45 


43.83 


<0.246 


0.988 


RQBL? 


A07 


104833.57+620305.0 


162.13990 


62.05139 




5.5 


19.85 




<0.281 


>1.330 


U? 


C05 


114748.99+351350.1 


176.95416 


35.23059 




2.5 


20.31 




<0.318 


0.833 


S? 


A07 


124225.38+642918.9 


190.60579 


64.48860 


0.042 


1.8 


17.12 


42.86 




>1.578 


G 


C05.VLA 


143139.15+600657.9 


217.91313 


60.11609 


0.413 


11.8 


20.21 


44.11 


<0.277 


1.246 


G 


A07 


150818.96+563611.2 


227.07903 


56.60314 


2.052? 


4.5 


19.51 


45.16 


<0.322 


>1.062 


U 


C05,VLA 


151115.49+563715.4 


227.81455 


56.62095 




0.6 


20.05 




<0.285 


>1.378 


Abs 


C05 


153304.11+453526.1 


233.26716 


45.59061 




3.7 


20.27 




<0.311 


1.215 


RQBL? 


A07.VLA 


154515.77+003235.2 


236.31574 


0.54311 


1.011? 


5.8 


18.81 


45.42 


<0.237 


>1.235 


U 


C05.VLA 


161004.03+253647.9 


242.51681 


25.61332 




6.5 


20.36 




<0.318 


0.876 


X 


A07.VLA 


165806.77+611858.9 


254.52822 


61.31638 


>1.410? 




20.69 


>46.41 


<0.025 


>1.747 


Abs 


C05 


211552.88+000115.4 


318.97035 


0.02097 




12.6 


19.47 






>1.192 


WLQ 


C05,VLA,Chandra,LC 


212019.13-075638.3 


320.07972 


-7.94400 




26.2 


19.86 




<0.242 


>1.168 


S? 


C05.VLA 


213950.32+104749.5 


324.95968 


10.79711 


0.296 


14.9 


20.11 


43.59 




0.985 


G 


C05 


224749.57+134248.1 


341.95658 


13.71338 


1.175? 


6.2 


18.27 


45.73 




>1.293 


WLQ 


C05.VLA, Chandra 


231000.83-000516.2 


347.50346 


-0.08785 




4.2 


19.00 




<0.192 


>1.254 


U? 


C05,VLA,LC 


232428.43+144324.3 


351.11848 


14.72344 


1.410 


1.2 


18.76 


45.65 




>0.844 


WLQ 


C05.VLA, Chandra 


a Tentative redshifts marked with ' 


?'. All objects have red 


shift upper limits z < 2.2, del 


ived by the fact that the Lyman forest is not seen in 


their SDSS spectra. If 


an object lacks a redshifi 


, we assume 
We considci 


z — 0.3 when 
proper motion 


calculating a^ rst and ck^ ss . 
measures to be significantly large 


f m > 30 n 




In milli-arcsec yr~ , 


fromlMunn et al.l d2004Ti. 


Lilli-arcscc yr because 



only 1% of spcctroscopically confirmed SDSS quasars have such large proper motion measures (sec Scction l4.2t . 

c Point spread function magnitude in the SDSS r filter. 

d At rest frame 5000 A. 

c Our classification (sec Sections^ and [5j. Radio-quiet weak-featured AGN arc classified as radio-quiet BL Lac objects (RQBL?), low-redshift 

WLQs (WLQ), or unknown (U). Contaminants (i.e., objects unlikely to be AGN with intrinsically weak lines) arc classified as absorbed AGN 

(Abs), galaxies (G), radio-loud BL Lac objects (RLBL), stars (S), or they arc unlikely the proper optical counterpart to a RASS X-ray source 

(X). Question marks denote low-confidcncc classifications. 

C05: listed a s a rad io-quiet BL Lac candidate in ICollinge et al.l ]2005l) ; A07: X-ray selected BL Lac candidate lacking radio emission in 
Anderson et al. (2007); VLA: we took follow-on 8.4 GHz radio observations with the VLA (see Section l3.1l >; Chandra: we took follow-on X-ray 
observations with Chandra (see Scction l3.2t ; LC: we have a light curve from SDSS Stripe 82 (sec Scction l3.3t . 

s pectrum into an ellipti cal galaxj[3 (using the template 
of lMannucci et al"]|2001l ) and a power law to meas ure each 
object's decomposed AGN optical flux (see lPlOf ) and its 
optical spectral index a (all local spectral indices are de- 
fined as /„ ~ v~ a "). We use these a measures, and we 
assume radio and X-ray spectral indices of a r = —0.27 
and a x = 1 .25, respectively, to calculate a ro and a ox at 
rest-frames 5 GHz, 5000 A, and f kcV. Limits on radio 
luminosities are estimated for sources not detected in the 
radio by FIRST/NVSS by assuming radio flux densities 
fiAGHz < 0.25 + 5a rms mJy, where 0.25 is the CLEAN 
bias and a rms is the uncertainty of t he FIRST survey a t 
the source's location on the sky (sec Bec ker et al.l fl995). 
We derive flux density limits from FIRST rather than 
NVSS because of its higher angular resolution. The op- 
tical magnitude in the SDSS filter closest to rest frame 
5000 A is used for calculating optic al lum inosities (and 
we correct for extinction with the Schlcgcl et al.l 119981 
dust maps). X-ray count rate limits are estimated for 
X-ray undetected sources as 6 counts divided by the ex- 
posure time of the nearest RASS X-ray sourceE3 We 



then calculate X-ray luminosities (or limits) at 1 keV 
using the Portable, In teractive Multi-Mission Simulator 
(PIMMS, Mu kalll993f) corr ecting for Galactic extinction 
with the lStarket al.1 (J1992T ) H I maps. 

Throughout, we consider objects with a ro < 0.2 to be 
radio-quiet for two reasons: 1) this is consistent with the 
traditional definition for radio-quiet AG N, that is, radio 
to optical flux rat ios R < 10 (Kellerm ann et al.l 119891 : 
iStockeet all 11992ft . 2) There are 285 optically selected 
SDSS BL Lac objects from lPlOl with both X-ray and ra- 
dio detections in RASS and FIRST/NVSS; these follow 
a Gaussian distribution in a ro with (a ro ) = 0.42 ± 0.08. 
Thus, any object with a ro < 0.2 is significantly differ- 
ent from a BL Lac object in its radio properties at the 
> 2.75cr level (>99.7% one-sided significance). This im- 
plicitly accounts for the non-simultineity of our radio and 
optical observations, and the fact that we assume the 
same radio spectral index for every object. We note this 
is slightly conservative because SDSS BL Lac objects de- 
tected by RASS are among the X-ray brightest and are 
radio- weaker on average (see lPlotkin et aLl l2008'). 



10 BL Lac obje cts are pro bably exclusively hosted by elliptical 
galaxies (e.g., sec Urrv et al. 2000). Host galaxy contamination 
can typically be seen if the host galaxy accounts for >20% of the 
observed flux. 

11 For the undetected RASS sources, we did not find adequate 



3. MULTIWAVELENGTH PROPERTIES 
3.1. Radio Observations with the VLA 



archived pointed X-ray observations (i.e., from SWIFT, Chandra, 
XMM-Newton, etc.). 
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We obtained radio observations with the VL A for 
17/26 radio-quiet BL Lac candidates from lC05 l and lAOfl 
These objects were targeted because they are among 
the optically brightest candidates (g < 21.3), and their 
FIRST/NVSS radio flux limits are not sensitive enough 
to place them in the radio-quiet regime (or they do not 
fall within the FIRST foo tprint). Due to scheduling con- 
straints, 5/26 IC05MA07I g < 21.3 radio-quiet BL Lac 
candidates falling between 07 fe < RA < 12^ were not 
observed. The remaining 4/26 IC05I/IA07I objects were ei- 
ther too optically faint (g > 21.3) to efficiently obtain 
improved a ro constraints with the VLA, or they were al- 
ready known to have a ro < 0.2 from their FIRST radio 
flux limits. 

VLA observations were taken in the D-array configura- 
tion on March 30, 2007 for five targets, and in the DnC- 
array configuration on June 8, 2008 for the other 12. All 
observations were performed in the X-band (8.4 GHz) 
with 26 antennas, with a total of 20-30 minutes integra- 
tion on each source (using 10 sec integrations). We used 
2 x 50 MHz intermediate frequencies (IFs) centered at 
v = 8435 MHz and 8485 MHz. Phase calibrators were ob- 
served every 10 minutes, and flux calibrations were tied 
to either 3C 48 or 3C 286. A summary of the observa- 
tions is given in Table[5J We observed at 8.4 GHz because 
many BL Lac objects are brighter at higher radio fre - 
quencies (i.e., a r ~ —0.27 typically. iStickel et al.lll99ll) . 
This combined with the VLA X-band's superior sensi- 
tivity compared to other observing bands makes seeking 
radio detections for BL Lac candidates most efficient at 
8.4 GHz. 

We reduced the data using standard routines in AIPS, 
flagging anomalous UV data when necessary. We cleaned 
the data and created maps with the IMAGR task, using 
a restoring beam ^8" and 2" x 2" pixels. We estimate 
background rms noises with TVSTAT, which range from 
&rms = 0.032 — 0.048 mJy beam -1 . We require radio de- 
tections to have 8.4 GHz flux densities > ha rms within 2" 
of the SDSS source position. We detect 4/17 sources in 
the radio, and we measure their flux densities with TV- 
STAT using apertures 8-10 pixels on a side. The largest 
separation between the radio and SDSS optical positions 
for these four radio detections is 1.4". 

We use the 8.4 GHz flux densities (or limits) combined 
with SDSS optical fluxes to determine a"j,°, which is 
again calculated between rest-frames 5 GHz and 5000 A. 
We follow the same methodology as in Section [2J placing 
5oy ms upper limits on the radio flux densities for the non- 
detections. A conservative 5<r rms threshold is used for 
consistency with upper limits placed by the FIRST sur- 
vey, and because an extraordinary claim of a radio-quiet 
BL Lac discovery requires very stringent constraints. We 
include the values (or limits) of the radio-optical broad- 
band spectral indices derived from our deeper VLA ob- 
servations (a"o°) m Table [2 and we also include the 
broad-band spectral indices derived from FIRST-SDSS 
{aH rat ) and RASS-SDSS (a™ ss ) for reference. 

We find only a single object (SDSS J0241+0009) to be 
radio-loud: it has an 8.4 GHz flux density of 2.44 mJy, 
and we calculate a^ a = 0.318. This object also appears 
in FIRST with a 1.4 GHz flux density of 0.67 mJy (below 
the nominal FIRST flux limit) and a£ rs *=0.249. This 
source did not appear in the FIRST catalog at the times 



IC05I and IA07I were published, which is why we included 
it in our VLA program. The discrepancy between a£* rs * 
and a"o° is primarily due to our choice of local radio spec- 
tral index a r = —0.27. Combining our 1.4 and 8.4 GHz 
flux density measurements, we find a r = —0.7. If we 
adopt a r = -0.7, then a^ st = 0.273 and a v r l a = 0.277. 
Regardless, SDSS J0241+0009 is definitively radio-loud. 
It also has a flat radio spectrum, it is detected in the 
X-rays by RASS, and we show in Section l3~3l that it is 
variable. Although it is on the relatively radio-weak tail 
of the radio-loud BL Lac distribution, its radio-loudness 
is not extraordinarily small. It should therefore be con- 
sidered a normal radio- loud BL Lac object. 

The other three VLA radio detections have a* 1 " val- 
ues that firmly place them in the radio-quiet regime. We 
note that one source (SDSS J2115+0001) might actu- 
ally be a blend of two radio point sources (it is not clear 
with the spatial resolution of our observations). Never- 
theless, we consider this source as a radio detection, and 
we report the total combined flux of the (potentially) 
two components. Even the combined radio flux confirms 
SDSS J2115+0001 is radio-quiet. 

We conclude from the radio observations that the ma- 
jority of weak-featured radio-quiet objects are indeed 
atypical in their radio properties from the much larger 
SDSS radio- loud BL Lac population. Of the 17 ob- 
jects with deeper VLA coverage, we only classify a single 
source as a confident radio-loud BL Lac object (SDSS 
J0241+0009). 

3.1.1. VLA Observations of High-Redshift WLQs 

Four z > 2.2 WLQs included in lC05F s list of poten- 
tial radio-quiet BL Lac candidates were included in our 
VLA programs, and none was detected in the radio. We 
list their 8.4 GHz radio properties in Table [3] for com- 
pleteness. Of th ese four objects, three have deeper X-ray 
observations in iShemmer et al.l (|2009T ). and we include 
those objects' a ox measures derived from their deeper 
X-ray observations J3 

3.2. X-ray Observations with Chandra 

Of the 17 radio-quiet BL Lac candidates with deeper 
VLA observations, 14 only have lower limits on a ox from 
RASS. Their RASS limits are not sensitive enough to 
determine if they are dissimilar to BL Lac objects in 
the X-ray, and deeper X-ray observations are required to 
compare their X-ray properties to the larger radio-loud 
BL Lac population. 

We observed three of these 14 objects in the X-ray with 
Chandra during Cycle-10. The Chandra targets were cho- 
sen as the optically brightest of the five sources with VLA 
observations taken in 2007 (which were the only VLA ob- 
servations we had in hand at that time). Chandra data 
were taken with the Advanced CCD Imaging Spectrom- 



12 Shcmmcr ct al. (2009) use slightly different parameters for 
calculating a ro and a ox . Among the differences, their broad-band 
spectral indices are defined with opposite sign compared to ours, 
and they are referenced at 2500 A in the optical and 2 keV in the X- 
ray. They also assume local spectral indices of 0.5 in both the radio 
and the optical, and they use a v = 1.0 in the X-ray. We convert 
their broad-band spectral indices, a! ro and a' ox , to be consistent 
with our adopted reference frequencies and local spectral indices 
using approximate transformations a ro = — 1.059oc' ro + 0.054 and 
a ox = -1.116a' - 0.058. 
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TABLE 2 
VLA Observations at 8.4 GHz of Radio-Quiet BL Lac Candidates 



SDSS Name 


Rcdshift 


Config a 


&rms 


fvla 


ccf* 


4T st 


a ox SS 


(J2000) 






(mjy beam - ) 


(mJy) 








012155.87-102037.1 


0.470 


DnC 


0.033 




<0.038 


<0.234 


>1.223 


013408.95+003102.4 




DnC 


0.047 




<0.123 


<0.266 


>1.152 


020137.70+002535.0 




DnC 


0.032 




<0.058 


<0.235 


>1.238 


024156.37+004351.6 


0.990 


DnC 


0.036 


0.43 


0.109 


<0.217 


>1.074 


024157.36+000944.1 


0.790? 


DnC 


0.041 


2.44 


0.318 


0.249 


0.863 


025046.47-005449.0 




DnC 


0.033 




<0.095 


<0.283 


>0.931 


025612.47-001057.8 




DnC 


0.032 




<0.111 


<0.288 


>0.876 


124225.38+642918.9 


0.042 


DnC 


0.035 




<-0.063 




>1.578 


150818.96+563611.2 


2.052? 


DnC 


0.032 




<0.128 


<0.322 


>1.062 


153304.11+453526.1 




DnC 


0.036 




<0.133 


<0.311 


1.215 


154515.77+003235.2 


1.011? 


DnC 


0.041 




<-0.009 


<0.237 


>1.235 


161004.03+253647.9 




DnC 


0.047 




<0.152 


<0.318 


0.876 


211552. 88+000115. 4 b 




D 


0.041 


0.32 


0.103 




>1.192 


212019.13-075638.3 




D 


0.041 




<0.065 


<0.242 


>1.168 


224749.57+134248.1 


1.175? 


D 


0.041 


0.98 


0.088 




>1.293 


231000.83-000516.2 




D 


0.040 




<0.035 


<0.192 


>1.254 


232428.43+144324.3 


1.410 


D 


0.048 




<0.011 




>0.844 



a Configuration of the VLA during observations. D-array data were taken on March 30, 2007 as VLA 
program AP524; DnC-array data were taken on June 8, 2008 as VLA program AP551. 
Might be a blend of 2 radio sources, both ~l" from the optical position. 



TABLE 3 
VLA Observations at 8.4 GHz of High-Redshift WLQs 



SDSS Name 
(J2000) 



Rcdshift Confie 



*rms _ fvla 

(mjy beam ) (mjy) 



031712.23-075850.3 2.699 

121221.56+534127.9 3.190 

123743.09+630144.9 3.535 

142505.59+035336.2 2.248? 



'nC 


0.035 


>nC 


0.038 


>nC 


0.044 


>nC 


0.038 



<-0.075 ■ ■ ■ 1.705° 

<-0.108 <0.072 2.095 b 

<-0.113 <0.058 >1.671 c 

<-0.058 <0.126 >1.236 d 



a Configuration of the VLA ckiring observations. Data were taken on June 8, 2008 as VLA program 
AP551. 

a ox value taken from Chandra detection in Shcmmcr ct al. (2009), with their reported a oa; value 

converted to conform to our a ox definition (sec footnotc ll2L 

c ol ox limit taken from XMM-Newton non-detection in Shcmmcr ct al. (2009). with their reported a ox 
limit converted to conform to our a ox definition (sec footnote ll2t . 

a ox limit taken from RASS non-dctcction. 



eter (ACIS; rGarmire et al.| [2003) at the nominal S3 aim- 
point, using faint telemetry mode. Exposure times of 
9.5 ksec, 3.3 ksec, and 5.0 ksec were achieved for SDSS 
J2115+0001, SDSS J2247+1342, and SDSS J2324+1443, 
respectively. The data were reduced using standard rou- 
tines in CIA0J3 The Chandra observations are summa- 
rized in Table |U 

We filter the X-ray observations from 0.5-6.0 keV to 
reduce the background. We then measure the number of 
X-ray counts within 6 pixel circular apertures centered on 
each source's position with the CIAO routine dmextract. 
The 6 pixel apertures include >90% of the encircled en- 
ergy at the S3 aimpointo The number of background 
counts are estimated over circular annuli with inner and 
outer radii of 10 and 20 pixels respectively. No photons 
are detected for SDSS J2115-t-0001; we place a n upper 
limit of 5 photons (99% confidence interval, see iGehreld 
119861) for SDSS J2115+0001, and we estimate a count 
rate < 0.0005 counts s" 1 . SDSS J2247+1342 is detected 



with 301^ photons (0.0091+^^5 counts s" 1 ); a weak 
detection is obtained for SDSS J2324+1443, with Q±\ 

photons (0.0012±g:ggg| counts s^H 

X-ray fluxes from 0.5-6.0 keV are estimated from the 
above count rates using PIMMS, correc ting for Galactic 
absorption with the IStark et al.l (|1992| ) H I maps, and 
assuming an X-ray spectral index a v = 1.25. We also 
estimate X-ray flux densities at 1 keV. Using the SDSS 
filter closest to rest-frame 5000 A, we find a^, andra > 



1.90, aZ andra = 1.55, and af x andra = 1.80 for SDSS 
J2115+0001, SDSS J2247+1342, and SDSS J2324+1443, 
respectively. 

The broad-band spectral indices derived from Chan- 
dra for these three objects place them on the very X- 
ray faint tail of the radio-loud BL Lac a ox distribution. 
Among the 637 radio- loud optically selected BL Lac ob- 
jects in lPlOl 294 are detected in the X-ray by RASS, with 
(oiox) = 1-11 ± 0.21 (note, this implicitly accounts for 
variability). The largest «™ ss for an object with a RASS 



13 Chandra Interactive Analysis of Observations, 
http : //exc . harvard . edu/ciao/ . 



Sec 



la See the Chandra P roposers ' 
http : //exc . harvard . edu/proposer /P0G/ 



Observatory Guide; 



15 Quoted err ors are 84.13% one-sided upper and lower confi- 
dence limits from Getircls (1986). These are calculated using Pois- 
son statistics, and they correspond to Gaussian ±lcr confidence 
intervals. 
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TABLE 4 
















Chandra 


BSERVATIONS 










SDSS Name 


Redshift 


Exp. Time 


N H a 


Count Rate 


fx c 


a aha ndra 


»& 


a o" SS 


(J2000) 




(kscc) 


(cm" 2 ) 


(counts s~ ) 


(erg s cm ) 








211552. 88+000115. 4 d 




9.5 


5.83 x 10 20 


<5xl0~ 4 


<3.2xl0~ 15 


>1.905 


0.103 


>1.192 


224749.57+134248.1° 


1.175? 


3.3 


4.58 x 10 20 


QltUxlQ- 4 


53.4xl0~ 15 


1.552 


0.088 


>1.293 


232428. 43+144324. 3 f 


1.410 


5.0 


4.15 x 10 20 


12+| xl0~ 4 


7.0xl0~ 15 


1.796 


<0.011 


>0.844 


a Column density from 


Stark et a] 


(1992) 















' From 0.5 - 6.0 keV, assuming a photon power law index T — 2.25. 

Corrected for Galactic absorption. 
1 Observation taken on Dec 24, 2008; Chandra ObsID 10388 
' Observation taken on Aug 7, 2009; Chandra ObsID 10387 

Observation taken on May 31, 2009; Chandra ObsID 10386 



X-ray detection is 1.66; only two BL Lac objects lacking 
RASS X-ray detections have larger a r % ss limits. If these 
three objects are typical BL Lac objects with abnor- 
mally weak radio-emission, then we might have expected 
their spectral energy distributions to peak at extremely 
high frequencies (i.e., they would be extreme high-energ y 
peaked BL Lac objects, see lPadovani fc Giommi|[l995af ). 
However, in this scenario we would expect them to be 
relatively X-ray bright, which is excluded by our Chan- 
dra observations. We thus conclude that the objects with 
deeper Chandra observations are dissimilar to BL Lac ob- 
jects in their X-ray properties, but the multiwavelength 
properties of these three objects are comparable to high- 
redshift SDSS WLQs (see Section Gp). 

Finally, we calculate Aa ra , the difference between the 
measured a^ andra values and that expected from their 
optical luminosities given the empirical a ox -L v (2500 A) 
relation in Equation 3 of I Just et all (|2007f ). L„(2500 A) 
is the specific luminosity at 2500 A rest-frame, and for 
our Aa ox measurements we re-calculate a^, andra at rest- 
frames 2 keV and 2500 A for consistency with lJust et al.l 
p00l . We find Aa ox =-0.052 (-0.26a) and 0.197 
(0.994cr) for SDSS 2247+1342 and SDSS J2324+1443, 
respectively. The values in parentheses are normal- 
ized to the standard deviation in a ox , a — 0.198, for 
AGN with optical luminosity dens ities between 10 30 - 
10 31 erg s- 1 Hz" 1 (see Table 5 of ISteffen et alJl2006h . 
The lack of a spectroscopic redshift (or an X-ray de- 
tection) for SDSS J2115+0001 precludes us from esti- 
mating its Att OI value. The measured Aa m for SDSS 
2247+1342 and SDSS J2324+1443 are not significantly 
different from zero, and we thus conclude that their X- 
ray properties are similar to normal radio-quiet quasars 
and WLQs (see lShemmer et al.ll2009h . 

3.3. Optical Flux Variability 

Finally, we consider flux variability, as we expect 
strong variability if these are in fact BL Lac objects. 
The SDSS has observed a ^300 deg 2 region of the sky 
(20 h 34 m < ajaooo.o < 04 h 00 m , |<J J20 oo.o| < 1-266°) in 
all five SDSS filters multiple times every Fal l since 1998. 
This region is referred to as Stripe 82 (see llvezic et al.l 
2007; Sesar et al. 2007). Seven of our low-redshift radio- 
quiet BL Lac candidates lie inside Stripe 82, each one 
with 50 - 70 epochs of observations (see Figure [1]) . The 
u and z filters tend to be noisier, so we do not consider 
them here. Most objects only have 1-2 epochs of observa- 
tions per year prior to Fall 2001, so we do not show these 



data in Figure [IJ but we do include pre-2001 epochs in 
the following analysis. 

We use the reduced Xr statistic to assess the level 
of flux variability in a given filter 03 If an object has 
Xr > 3, then we co nsider it to show variability (e.g., 
see iSesar et al.l l2007f) . We find five of the seven sources 
to be variable in the g filter (SDSS J0241+0043, SDSS 
J0241+0009, SDSS J0256-0010, SDSS J2115+0001, and 
SDSS J2310-0005). These five objects show flux varia- 
tions over short and long timescales (ranging from days 
to years), and flux variations tend to be correlated be- 
tween filters. We note that one of the variable objects 
is SDSS J0241+0009, which we classify as a radio-loud 
BL Lac object based on its radio and X-ray properties. 

The measured values of Xr an d the rms scatter, a rms , 
for each source in each filter are listed in Table [SJ We 
include a^ a and a™ ss for reference, and we list the aver- 
age measured photometric error for each source in each 
filter, (a m ). Observations typically have photometric un- 
certainties of ^0.03 mag, and the largest degree of vari- 
ability we detect is a rms ~ 0.2 mag (SDSS J0256-0010). 
The two objects that do not appear variable are likely 
stars. Neither has a spectroscopic redshift, and one of 
them (SDSS J0134+0031) has a very large proper mo- 
tion measure ([J, ~ 150 milli-arcsec yr ). 

The detection of variability does not necessarily indi- 
cate that these are BL Lac objects, because radio-quiet 
quasars are also variable. We now compare the level 
of variability of the radio-quiet BL Lac candidates to 
normal radio-loud BL Lac objects and to radio-quiet 
quasars. There are 14 optically-selected BL Lac can- 
didates in P10 with Stripe 82 light curves that are de- 
tected in the radio by FIRST/NVSS and are radio-loud; 
we calculate xf f° r those objects in the g-filterQ We 
also examine a sample of ^2000 spectroscopically con- 
firmed radio-qu i et SPS S DR5 quasars in Stripe 82 from 
ISchneider et al.l (J2007I ). and we measure xl f° r eacn of 
these radio-quiet quasars. These ^2000 quasars are se- 
lected because: 1) they have Stripe 82 light curves with 



16 2 

A. r 



Ef [(mi 



(m))/o"i] /(N — 1), where m; and en are 
the observed magnitudes and uncertainties, respectively, in a given 
filter, (m) is the average magnitude, and N is the number of epochs 
in that filter. 

17 Because SDSS J0241+0009 has radio-loud a ro measures from 
both FIRST and our follow-up VLA observations, we do not con- 
sider it a radio-quiet BL Lac candidate in th e foll owing. It was 
recovered as a radio-loud BL Lac candidate in lPldl . and it is thus 
included in the 14 object subset of radio-loud BL Lac objects with 
Stripe 82 coverage. 
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Fig. 1. — Stripe 82 light curves in the g, r, and i filters for seven radio-quiet BL Lac candidates; data taken prior to Fall 2001 are not 
shown. The solid lines illustrate the average magnitude of each source, and the dotted lines mark ±3 times the average photometric error. 
All panels span 1 mag along the y-axis. The two non-variable objects (with \r < 3) are likely stars. 



at least 10 epochs of data in each filter; 2) they are in 
the FIRST survey's footprint and have a ro < 0.2, de- 
rived from FIRST radio detections or limits; (3) they 
have g > 18.7, so that their optical fluxes are similar to 
our radio-quiet BL Lac candidates with Stripe 82 cov- 
erge; and (4) they have 0.5 < z < 2.2. This rcdshift 
range encompasses the redshifts of our seven radio-quiet 
BL Lac candidates in Stripe 82. All have z < 2.2, since 
we do not see the Lya forest in their SDSS spectra; those 
lacking spectroscopic redshifts probably have z > 0.5, or 
else we would likely see host galaxy spectral features in 
their optical spectra fsee lPlOf ). We restrict our compari- 
son quasar sample to these redshifts so that our Xr mea- 
sures are calculated at approximately similar rest-frames. 
For simplicity, the following analysis is illustrated with 
values derived in the g filter, but our conclusions do not 
change if we use a different filter (or if we choose filters 
to cover similar rest-frame wavelengths). 

Figure[2]shows Xr vs - a ro for the six radio-quiet BL Lac 
candidates with Stripe 82 light curves (squares, not in- 
cluding SDSS J0241+0009) and for the 14 radio-loud 
BL Lac candidates (circles) from lPldl that lie in Stripe 82. 
We use VLA data when available to calculate a ro . The 
^2000 radio-quiet quasars have an average (Xr) = 65, 



and 84% (i.e., the one-sided 1<t confidence interval) have 
Xr < 106, and 97.5% (i.e., the one-sided 2ct confidence 
interval) have Xr < 245 (marked by the dotted lines in 
Figure [2]) There is a marked deficit of extremely vari- 
able radio-quiet objects. 

The level of flux variability is likely dominated by dif- 
ferent physical processes in radio-quiet and radio-loud 
quasars. Using a sample of ^100 quasars covering a wide 
range of redshifts, optica l luminosities, and radio loud- 
nesses, iKellv et all (|2009l ) suggest the flux variability of 
radio-quiet quasars is cause d by a stochastic proce ss; this 
conclusion is confirmed by Ma cLeod et al.l (l2010f). O ne 
potential mechanism proposed by IKellv et al.l (|2009[ ) to 
explain the observed variability for radio-quiet quasars is 
turbulent magnetic fields within the accretion disk. For 
radio-loud quasars, however, additional mechanisms re- 
lated to the radio jet need to be invoked. Here, variabil- 
ity is commonly modeled as sho cks propagating through 
an inhomogeneous iet (e.g.. see iMarscher fc Gearl 119851 : 
lUlrich et aJJll997l: lAUeret al.lll999l) . ~ 

The Xr values in Figure suggest that the mechanism 
driving variability in the radio-quiet BL Lac candidates 



in 



Only 19 quasars have \r < 3 in the g filter. 
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TABLE 5 
Optical Variability of Seven Objects in Stripe 82 



SDSS Name Redshift 

(J2000) 



l ro ^oa; 
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(mag) 



(Om) 

(mag) 
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<0.123 


>1.152 
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0.033 
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1.221 


0.032 


0.028 
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0.777 


0.029 


0.034 
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0. 
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.074 
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r 
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9 
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<0. 


095 


>0 
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9 
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.905 a 


9 

V 

i 


13.577 
4.080 
2.988 


0.085 
0.040 
0.037 


0.025 
0.022 
0.022 


231000.83-000516.2 


<0. 


035 


>1 


.254 


9 
r 

i 


15.379 
19.540 
12.967 


0.050 
0.066 
0.062 


0.016 
0.015 
0.017 



X-ray flux limit from Chandra 
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Fig. 2. — Reduced \r °f A ux variability in the q filter vs. a ro . 
Circles: 14 radio-loud BL Lac objects from IP1CI with Stripe 82 
light curves. Squares: six radio-quiet BL Lac candidates with 
Stripe 82 light curves and a ro measures (or limits) smaller than 
0.2. The vertical solid line shows the traditional boundary between 
radio-loud and radio-quiet AGN. The horizontal dashed-dotted line 
shows Xr = 3, above which we conclude a source is variable. The 
dotted lines mark the maximum values of Xr we expect 84% (lcr) 
and 97.5% (2<r) of radio-quiet quasars to show. The variability of 
the radio-quiet BL Lac candidates is consistent with that expected 
for radio-quiet quasars, although there is some overlap between the 
high and low variability tails of the radio-quiet quasar and BL Lac 
distributions. This figure is similar if we rather plotted \r f° r the 
r or i filters. 



is similar to that for radio-quiet quasars. To test this as- 
sertion, we assume the observed distribution of Xr f° r the 
four variable radio-quiet BL Lac candidates is randomly 
drawn from the \r distribution of the ^2000 DR5 radio- 
quiet quasars. We ignore the two non- variable radio- 
quiet BL Lac candidates because they are likely stel- 
lar. We also assume the 14 radio-loud BL Lac objects' 
Xr measurements are randomly drawn from the DR5 
quasars. The above null hypotheses are then tested via 
Monte Carlo simulations, as described below. This test is 
designed to achieve the same purpose as a Kolmogorov- 
Smirnov test, but we use Monte Carlo simulations be- 
cause of the small sample sizes. 

Of the 14 IP10I radio-loud BL Lac objects with Stripe 82 
light curves, 12/14 have Xr > 65 (the mean of the 
DR5 quasars), 9/14 have \l > 106 ( lcr )> and 9 / 14 have 
Xr > 245 (2cr). To test the likelihood that the ~2000 
DR5 quasars are the parent population of the 14 radio- 
loud BL Lac objects, we randomly choose 14 Xr values 
from the DR5 quasar distribution. We then count the 
number of randomly selected Xr values larger than 65, 
106, and 245; we repeat 10 6 times. Our Monte Carlo 
simulations show there are p=0. 000091, 0.000061, and 
< 0.000001 chances of randomly drawing 12, 9, and 9 ob- 
jects with Xr values larger than 65, 106, and 245, respec- 
tively. We thus conclude the radio-loud BL Lac objects 
are statistically different from the DR5 quasars in their 
variability properties. 

For the four radio-quiet BL Lac candidates, no object 
has measured Xr > 65. We perform a similar test, this 
time randomly drawing only four Xr values from the dis- 
tribution of DR5 quasars, and we again repeat 10 6 times. 
We find there are p=0.19, 0.50, and 0.90 chances of ran- 
domly choosing zero quasars with Xr >65, 106, and 245, 
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respectively. Thus, we conclude it is plausible the level 
of variability displayed by our radio-quiet BL Lac candi- 
dates is similar to that of normal radio-quiet quasars. 

We of course caution that one should not draw 
strong statistical conclusions from these small sample 
sizes; however, our Monte Carlo simulations suggest 
the radio-quiet BL Lac candidates are similar to radio- 
quiet quasars, but the radio-loud BL Lac objects are 
not. This conclusion is strengthened after also con- 
sidering their radio and X-ray properties. It is also 
important to mention that there is overlap in \r be- 
tween the most variable radio-quiet quasars and the 
least variable radio-loud BL Lac objects. Thus, we 
cannot exclude that some of the more variable radio- 
quiet BL Lac candidates (e.g., SDSS J0256-0010) have 
weakly boosted continua, but that does not appear to be 
the dominant source of radiation. A similar conclusion 
was reached for some relatively radio-bright high-redshift 
WLQs (e.g., SDSS J1408+0 205 and SDSS J1442+0110; 
iDiamond-Stanic et al.l 12009). The most variable radio- 
quiet BL Lac candidates should be placed at high priority 
for polarimetric follow-up. 

4. RECOGNITION OF SOURCES UNLIKELY TO BE 
LINELESS AGN 

Equipped with the above observations, we now iden- 
tify objects that are unlikely to be AGN with intrinsi- 
cally weak emission lines. Here we discuss all 26 radio- 
quiet BL Lac candidates in Table HJ not just the 17 with 
deeper VLA observations. Improvements to the SDSS 
data reduction pipelines are implemented prior to each 
data rel ease, so we also inspected the DR7 spectra of the 
IC05I and lA07l radio-quiet BL Lac candidates (which were 
selected from earlier SDSS data releases). 

4.1. Normal Radio- Loud BL Lac Object 

SDSS J0241+0009 is a normal radio loud BL Lac ob- 
ject based on its radio and X-ray properties, and its op- 
tical variability. 

4.2. Stars 

Three of the 26 radio-quiet BL Lac candidates show 
large proper motion in Table [TJ which we define as 
fi > 30 milli-arcsec yr _1 (less than 1% of spectroscop- 
ically confirmed SDSS quasars have pro per motion mea- 
sures larger than 30 milli-arcsec yr _1 ; iSchneider et al.l 
2007). One of these measurements is spurious: SDSS 
J0241+0009 (fi — 54 milli-arcsec yr _1 ) is a very con- 
vincing radio- loud BL Lac object (see Section |4~T|) . The 
other two sources with large proper motion measures 
(SDSS J0134+0031 and SDSS J0201+0025) both have 
/i > 100 milli-arcsec yr _1 , they lack spectroscopic red- 
shifts, and they have point source optical morphologies [_j 
Neither is detected by our deeper radio observations or 
by RASS in the X-ray. We have an optical light curve 
for SDSS J0134+0031 (see Section O, and it does not 
show variability over an entire decade. Those two proper 
motion measures are likely not spurious, and those two 
objects are removed from the sample as potential stars. 
We also remove SDSS J0250-0054 as a likely star be- 
cause it does not show flux variability in its Stripe 82 

19 Both objects are from C05, and their proper motions were 
less certain at the time of publication. 



light curve, it does not have a measured redshift, and it 
is not detected in our 8.4 GHz VLA observations. The 
spectrum of SDSS J1147+3513 might be contaminated 
by a bright nearby star, and SDSS J2120-0756 might 
also show weak stellar absorption features in its DR7 
spectrum. We thus remove those two sources as well 
(neither has a measured redshift). In all, five of the 26 
radio-quiet BL Lac candidates are removed because they 
are probably stars. 

4.3. Galaxies 

Examination of the improved SDSS spectral reduc- 
tions of SDSS J1431+6006 and SDSS J2139+1047 show 
they have Ca II H/K breaks larger than 40%, and we 
re-classify them as elliptical galaxies. We also remove 
SDSS J1242+6429, which is probably a galaxy: its very 
low-redshift (z = 0.042) and low optical luminosity 
(\ogv L Vo < 43 erg s _1 ) suggest it is not an AGN. 

4.4. Absorbed AGN 

The weak emission lines from SDSS J1511+5637 and 
SDSS J1658+6118 might be explained by absorption. 
Both objects have extremely red spectra (optical spec- 
tral indices a u — 4.4 and 5.9, respectively), and SDSS 
J1658+6118 is ad ditionally classified a s a broad absorp- 
tion line quasar in lTrump et al.l (J2006I ). 

4.5. Misidentified Optical Counterparts to RASS X-ray 

Sources 

There are sometimes multiple optical sources within 
each 60" RASS X-ray circle, and we believe two radio- 
quiet BL Lac candidates (SDSS 0755+3525 and SDSS 
J1610+2536) are not the proper optical counterparts. 
SDSS 0755+3525 is 48" from the center of the RASS 
error circle. There is a more likely o ptical counterpart 
(identified as a starburst galaxy by iBrinkmann et all 
I2000D just 3" away from the RASS source and 0.2" from 
a FIRST radio source (FIRST J075551. 3+352635). For 
SDSS J1610+2536, there is a bright star 16" away (cat- 
aloged as V1024 Her, an eclipsing variable star with 
V = 12.5 mag in SIMBAD) that is also inside the RASS 
X-ray circle. V1024 Her is probably the X-ray emitter. 
Although its spectral type is not available in the liter- 
ature (and it is too bright to have an SDSS spectrum), 
its optical magnitude and RASS X-ray count rate are 
consistent wit h those exp ected from a K or M star (see 
Figure 3 of A giieros et al.ll2009l ) . We remove both sources 
from our list of radio-quiet BL Lac candidates. 

5. DISCUSSION 

We are left with 13 AGN with weak emission features 
(see Figure[3|), 10 of which are confirmed to be radio-quiet 
from our VLA observations. The other three sources' 
FIRST radio limits are not sensitive enough to determine 
if they are radio-quiet. Our 13 surviving weak-featured 
AGN are likely a heterogeneous combination of at least 
two populations of objects. Here, we classify them as 
likely low-redshift WLQs, as radio-faint BL Lac objects, 
or as uncertain. These classifications are listed in Ta- 
ble Q] We also compare them to other low- and high- 
redshift radio-quiet AGN in the literature with weak or 
absent spectral features. 
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Fi g. 3, — SPSS spect r a of th e 13 surviving weak-featured AGN in the observed frame. In the bottom-right panel we show, for reference, 
the Vanden Berk ct al. (2001) composite SDSS quasar spectral template (shown in the rest-frame), with prominent broad emission lines 
labeled. The expected locations of these lines are drawn for the radio-quiet BL Lac candidates with redshifts. Redshifts smaller than 
0.5 are derived from host galaxy absorption features (not labeled). Our classification as Weak Line Quasar (WLQ), Radio-quiet BL Lac 
Candidate (RQBL?), or Unknown (U) from Table[JJis given after each object's redshift. Some panels show a residual sky line at 5500 A. 



5.1. Low-redshift WLQs 

The multiwavelength colors of our 13 radio-quiet 
BL Lac candidates are shown in Figure 2J The three ob- 
jects with Chandra coverage (SDSS J2115+0001, SDSS 
J2247+1342, and SDSS J2324+1443, open upside down 
blue triangles) are radio and X-ray fainter than normal 
BL Lac objects. SDSS J2115+0001 does not have a red- 
shift, but it does have an 8.4 GHz radio detection im- 
plying a ro = 0.1. It also shows weak variability in its 



Stripe 82 light curve. Smi th et all (|2007f l perfor med a 
polarimetric survey of 42 BL Lac candidates from IC051 
and they did not detect strong polarization from either 
SDSS J2247+1342 or SDSS J2324+1443 (both have po- 
larization P < 1%). This further suggests they should 
not be unified with BL Lac objects. Based on their low 
polarization and multiwavelength colors, they may in- 
stead be low-redshift analogs to WLQs. 
The prototype high-redshift WLQ was discovered by 
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Fig. 4. — Broad-band spectral indices of 13 surviving z < 2.2 radio-quiet BL Lac candidates, with arrows denoting limits. Filled blue 
squares: two objects with RASS X-ray detections retained (with low-confidence) as potential very radio-faint BL Lac objects. Open 
blue diamonds: eight objects lacking sensitive X-ray constraints. Open upside-down blue triangles: three objects with Chandra X-ray 
observations, these are likely low-redshift WL Qs. For ref e rence , we show 14 high-redshift (2.7 < z < 5.9) WLQs (filled red triangles) 
with Chandra/XMM X-ray follow-up from Shemmcr et al. (2006) and Shcmmer ct al. (2009). Their radio information is typically derived 
from FIRST/NVSS, except for three WLQs with follow-up VLA observations in Table |3] We also plot ~600 optically selected radio-loud 
BL Lac objects. Circles denote BL Lac objects with an X-ray detection in RASS, and arrows mark X-ray limits for objects lacking RASS 
detections. The dashed line shows the traditional boundary between radio-loud/quiet quasars (a ro = 0.2). 

the detection of a relativistically boosted core for PG 
1407+265, which might (help) explain its odd spectral 
characteristics . Even with the detection of a radio core, 
iBlundell et all (|200l do not believe PG 1407+265 has 
b een radio loud within t he last 10 6 -10 7 years. 

iLeighly et al.l (|2007a[ ) present detailed multiwave- 
length spectral coverage of PHL 1811 (z=0.192), an un- 
usual X-ray weak narrow-line type 1 quasar. They detect 
relatively strong Ha and H/3 emission, but forbidden, 
semiforbidden, and high ionization lines are extraordi- 
narily weak or absent. PHL 1811 also shows relatively 
strong Fe II and Fe III emission. We note in particular 
that the SDSS spectra of SDSS J2247+1342 and SDSS 
J2324+1442 look similar to PHL 1811, and we might ex- 
p ect them to al s o have strong Fe II and Fe III lines. 

ILeighly et al.l (|2007al) run photoionziation models and 
find almost all of the odd spectral characteristics of PHL 
1811 can be explained by its soft SED, as fewer high- 
ionization species are formed than would result from a 
harder (and more normal) SED. Most of the radiative 
cooling is performed via Hydrogen lines (which can be ex- 
cited by the soft contin uum), hence the relat ively strong 
low-order Balmer lines. Leighly ct al. (2007b) show PHL 
1811 has a steep X-ray photon spectral index (r ~ 2.3 
over 0.3-5.0 keV), and they attribute this to an excep- 
tiona lly high (p erhaps super-Eddington) accretion rate. 
IShemmer et al . (2009) note that Chandra detected high- 
redshift SDSS WLQs do not appear to have unusually 
steep hard X-ray spectra, but they stress that X-ray ob- 
servations with improved photon statistics are required 
to test this properly. Infrared spectroscopy to determine 



IFan et al.l (fl999l SDSS J1532-0039; z = 4.62), who 
remarked on a "high-redshift quasar without emission 
lines." SDSS J1532-0039 was not det ected in the X- 
ray w ith follow-on Chandra observations (Shc mmer et al.l 
2006); it also remained undetected by the VLA (clearly 
placing it as radio-quiet), and optical observations 
also did not find any polarizati on or strong variabil - 
ity (|Diamond-Stanic et al. 2009). Anderson et al.1 1)20011 ) 
discovered two more z = 4.5 — 4.6 W LQs, and a fourt h 
discovery at z = 5.9 was rep orted by IFan et al.l (|2006| ). 
iDiamond-Stanic et al.l (2009) present a sample of ~70 
SDSS WLQs (z > 3), defined by Lya+N V REW < 
15 A. They show that WLQs constitute a prominent ex- 
cess of objects in the low Lyct+N V equivalent width tail 
of the high-redshift quasar distribution (which otherwise 
follows a log-normal distribution) , and there is no corre- 
sponding excess of objects in the high- REW tail. 

While the SDSS has produced a relatively large num- 
ber of high-redshift discoveries, there are only a hand- 
ful of low-redshift candidates in the lit erature. To our 
know ledge, PG 1407+265 (z = 0.94; iMcDowell et~aTl 
I1995D was the first AGN with weak emission features 
not interpreted as a BL Lac phenomenon. PG 1407+265 
has radio, X-ray, and o ptical flux ratios typica l of nor- 



1995) present 



mal radio-quiet quasars. IMcDowell et al.1 
UV, optical, and near-infrared spectra, and they find all 
lines (in particular, Lya, C IV and Mg II) have unusu- 
ally small equivalent widths. However, they do detect 
comparatively strong Ha by this odd source's standards 
(REW = 126 A, which is still relative l y wea k com- 
pared to most quasars). IBlundell et al.l (|2003|) report 
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the strength of Ha and/or H/3 for SDSS WLQs would 
also be illuminating, and H/3 li ne widths may allow ac- 
cretion rate estimates (e.g., see lShemmer et al"1l2004f) . 

A proposed lower-redshift analog to high-redshift 
WLQs was serendipitou sly discovered in t h e SPS S spec- 
troscopic database by iHryniewicz et all (|2010L SDSS 
J094533.99+100950.1; z = 1.66). SDSS J0945+1009 
has a normal optical continuum compared to radio- 
quiet quasars, but it shows very weak Mg II emission 
(REW ~ 15 A) and absent C III] and C IV emission. 
SDSS J0945+1009 is similar to the intermediate redshift 
radio-quiet BL Lac candidates presented in Table [T] how- 
ever, even its weak Mg II is too large to make it into our 
SDSS BL Lac samples. 

5.2. Radio- Quiet BL Lac Candidates 

Two objects (SDSS J1008+1149 and SDSS 
J1533+4535) are detected in the X-ray by RASS, 
with X-ray brightnesses typical of BL Lac objects 
(a ox =0.99 and 1.22, respectively, see Figure H]). Both 
objects appear to be the best optical counterparts in 
each RASS error circle, but of course Chandra observa- 
tions would be helpful for confirmation. These sources 
are very interesting, but we cannot currently assert they 
are extraordinary radio-quiet BL Lac objects: we do not 
have deeper VLA constraints for SDSS 1008+1149, and 
its FIRST flux limit (a ro < 0.25) does not firmly identify 
it as radio-quiet. Its redshift (z = 0.26) is also very 
typical of the parent SDSS BL Lac sample, and, similar 
to SDSS J0241+0009, this could simply be a relatively 
radio-faint but otherwise normal high-energy peaked 
BL Lac object. SDSS J1533+4535 was not detected in 
deeper VLA observations, and it is radio-fainter than 
normal BL Lac objects at the >3.4cr level. However, its 
SDSS spectrum is relatively noisy. Both sources should 
be placed at high-priority for further study. 

5.3. Unknown 

The remaining eight objects are very likely AGN, but 
their RASS X-ray limits are not sensitive enough to de- 
termine if they are X-ray weaker than BL Lac objects 
(see Figure HJ . Deeper X-ray observations are necessary 
to determine if they are best unified with BL Lac ob- 
jects or with WLQs. We only lack deeper VLA obser- 
vations for two objects - SDSS J0901+0314 and SDSS 
J1048+6203. 

Three objects (SDSS J0241+0043, SDSS J1508+5636, 
and SDSS J1545+0032) are particularly difficult to unify 
with BL Lac objects because they are at relatively high 
redshifts: they all have z > 1, and the me dian red- 
shift of radio-loud SDSS BL Lac objects from IPIOI with 
reliable re dshifts is 0.34. Polarimetric observations by 
ISmith et all ((20071) also found SDSS J1545+0032 to not 
be highly polarized (P ~1%). All three objects have 
df" < 0.2, and we suggest they are most likely lower- 
redshift WLQs (although we do so with low-confidence 
given their poor X-ray constraints). 

Two objects, SDSS J0121-1020 and SDSS 
J0901+0314, have redshifts (z = 0.47 and z = 0.46, 
respectively) typical of normal BL Lac objects. SDSS 
J0121-1020 is definitively radio-quiet (a ro < 0.04, 
>4.75u radio fainter than normal BL Lac ob- 
jects); we do not have deeper VLA observations 



for SDSS J0901+0314, but its FIRST radio flux limit 
(a(l rst < 0.20) indicates it is radio- fainter than BL Lac 
objects at the >2.75<r level. Deeper X-ray observations 
should be taken for both sources to determine if they 
are more likely very low-redshift WLQs or if they are 
very radio-faint BL Lac objects. 

Neither source shows obvious signatures of star forma- 
tion or for obscured emission lines, and their Ca II H/K 
depressions are much smaller than expected for normal 
elliptical galaxies. We do note that weak stellar absorp- 
tion lines are seen in the spectrum of SDSS J0901+0314, 
which suggests it might be a post-starburst galaxy. How- 
ever, an enhanced blue continuum from recent star for- 
mation alone is unlikely strong enough to explain its weak 
Ca II H/K break and lack of strong emission lines. For 
example, SDSS J0901+0314 shows only very weak US 
with REW of 1-2 A in absorption, and post-starburst 
galaxies typically require H<5 stronger than 5 A in ab- 
sorption (e.g., see Goto 200j|). 

Interestingly, these two objects appear remark- 
ably similar to the radio-quiet BL Lac candidate 
2QZ J215454 . 3-30 5654 (z = 0.49) discovered by 
ILondish et all (|2051 : 2QZ J2154-3056 shows only weak 
[O III] and perhaps some very weak Balmer emissionPI 
Both SDSS J0121-1020 and SDSS J0901+0314 also 
show weak forbidden lines from Oxygen and very weak 
narrow H/3 (all lines have REW < 5 A). Londis h et al.l 
(2004) do not detect optical polarization or radio emis- 
sion (a ro < 0.08), arguing against a beamed synchrotron 
continuum for 2QZ J2154— 3056. X-ray emission is not 
seen in RASS (fx < 10~ 14 erg s _1 cm -2 ), and we esti- 
mate from their Figure 2 that a ox > 1.2. They do detect 
marginal variability; the variability, however, is seen in 
only one of four epochs over a ~2-3 month period with 
a reduced xl = 4- The lack of relatively strong H/3, and 
its redder spectrum, makes 2QZJ 2154—3056 dissimilar 
to PHL 1811 and PG 1407+265. 

2QZJ 2154—3056 and the similar SDSS sources are un- 
likely explained by near- or super-Eddington accretion 
rates (as in PHL 1811), because then we would expect 
relatively strong H/3. It is alternatively possible to cre- 
ate lineless AGN with exceptionally low accretion rates 

(<10 -2 -10~ 3 MEdd), as has been invoked for so called 
"naked AGN" (i.e., Type 1 Se yfert galaxies lacking broad 
emission lines; Hawkins 2004) and X- ray bright, opticall y 
normal galaxies (XBONGS; e.g., see lTrump et alJl2009f l. 
At these low accretion rates, the inner region of the ac- 
cretion disk becomes radiatively inefficient. Thus, a radi- 
ati on pressure driven wind never forms th e BELR (e.g., 
see iNicastrol |2000l: iNicastro et al.l 120011 ). lElitzur fc Hoi 
(2009) reach a similar conclusion based on mass conser- 
vation arguments. 

SDSS J0121-1020 and SDSS J0901+0314 are, how- 
ever, unlikely to have low-enough accretion rates for their 
BELRs to vanish. Eve n if they have extre mely massive 
black holes (10 10 M L lElitzur fc Hoi (|2009ft predict they 
would require radiatively inefficient accretion flows and 
bolometric luminosities L\> i < 5 x 10 41 erg s _1 for their 
BELRs to disappear; this scenario is excluded by their 
optical luminosites. Also, in the absence of a radiatively 

20 Ha is not in the spectral coverage of their optical spectrum. 
Ca II H/K and G-band absorption from the host galaxy are also 
detected. 
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driven wind, the outflow is likely replaced by a jet, so 
we would expect very low accretion rate syste ms to also 
have strong radio emission (e.g., see iHol 120081) . 

It would be surprising if SDSS J0121-1020, SDSS 
J0901+0314, and 2QZJ 2154-3056 were relatively radio- 
faint but otherwise normal BL Lac objects with their 
emission lines simply d i luted by starlight from the host 
galaxy. iMarcha et al.l (|1996l ) argue that the strongest 
emission features in weakly beamed BL Lac objects show 
equivalent widths up to ^50-60 A when the equiva- 
lent width measures are referenced to the contaminat- 
ing starlight continuum. In those cases, a boosted syn- 
chrotron continuum still acts to weaken the lines slightly 
(albeit less than if the AGN was more highly beamed). 
We do not believe the host galaxy on its own could ren- 
der Ha or H/3 completely undetectable. In the absence 
of a boosted continuum, Ha and H/3 in these radio-quiet 
BL Lac candidates are likely intrinsically weak, but the 
mech anism is currently unclear (but see lLondish et al.1 
120041 for potential explanations). 

Finally, three other objects (SDSS J0256-0010, SDSS 
J1048+6203, and SDSS J2310-0005) remain as weak- 
featured AGN, but we retain them with lower-confidence 
because they have relatively noisy SDSS spectra. All 
three also lack spectroscopic redshifts, but SDSS 
J0256-0010 and SDSS J2310-0005 show variability in 
their Stripe 82 light curves. Both of these sources are 
also confirmed to be radio-quiet in their deeper VLA 
observations. We do not have VLA data for SDSS 
J1048+6203, and its FIRST limit is not terribly sensi- 
tive (a ro < 0.281). It is possible this source is a normal 
radio-loud BL Lac object, but deeper radio and X-ray 
observations are necessary to settle this issue. 

5.4. Radio- Quiet BL Lac Candidates from lPlul 

There are 71 additional z < 2.2 o bject s identified 
as radio-quiet BL Lac candidates in IPIOJ All have 
FIRST/NVSS radio flux densities or limits in the radio- 
quiet regime (a ro < 0.2), so they are not just opti- 
cally faint BL Lac objects with radio fluxes below the 
FIRST/NVSS limits. All but five of their objects have 
extragalactic redshifts, so th ere appears to be little con- 
tamination from stars. IPIOI argue many of these objects 
are distinct from BL Lac objects because they tend to 
be at systematically higher redshifts and are more lumi- 
nous than radio-loud BL Lac objects. However, deeper 
X-ray observations are necessary: the majority (95%) 
do not have RASS X-ray detections, and their RASS 
flux limits are not sensitive enough to discriminate them 
from BL Lac objects in the X-ra y. Based on our con- 
clusions from our analysis of the IC05I/IA07I sample, we 
ex pect deeper X-ray observations will show the majority 
of P10 objects to be low-redshift WLQs. 

6. CONCLUSIONS 

Here, we present follow-up radio and X-ray observa- 
tions for a subset of z < 2.2 radio-quiet BL Lac candi- 
dates from a sam ple of 26 sources originally presented 
in IC05I and IA07I . When available, we also consider flux 
variability and proper motion measures, and we conclude 
that 13 of the 26 objects are AGN with weak-emission 
lines. All but three of these 13 are definitively radio- 
fainter than BL Lac objects at the >2.75cr level (the 
other three objects are not detected in the radio, but 



their radio flux limits are not sensitive enough to de- 
clare them radio-quiet), making it difficult to interpret 
them as beamed radio galaxies. We thus confirm that 
the SDSS has recovered a population of radio-quiet AGN 
with intrinsically weak emission lines. That is, relativis- 
tic beaming does not appear to be the dominant cause 
of these objects' featureless spectra (but we also cannot 
exclude the possibility that weak jets exist in some of 
these sources). These 13 weak- featured AGN are likely a 
heterogenous mix of at least two populations of objects. 

Of the 13 AGN, three with deeper Chandra X-ray 
observations are likely low-redshift WLQs, analogs to 
the ^10 2 object high-redshift (z > 2.2) population dis- 
covered by the SDSS. Other low-redshift WLQ candi- 
dates have_j2e^n_djscoyered serendipitojjsly in th e SDSS 
(e.g.. lHrvniewicz et alJ "2010). but the IC05MA071 BL Lac 
surveys were the first systematic searches of the SDSS 
database capable of recovering low-redshift WLQs. 

There are two radio-quiet BL Lac candidates that have 
X-ray detections in RASS and are too X-ray bright to be 
unified with WLQs, but their a™ ss measures are similar 
to those of BL Lac objects. We retain these as radio- 
quiet BL Lac candidates, but with low-confidence. One 
object's radio flux limit is not sensitive enough to declare 
it radio-quiet (i.e., it could be a relatively radio- faint but 
still normal radio- loud BL Lac object). The other ob- 
ject is definitively radio-quiet, but its optical spectrum 
has relatively low S/N. The fact that we do not find 
a single high-confidence radio-quiet BL Lac object sup- 
orts the notion th at all BL Lac objects are radio-loud 
Stocke et al1ll990h . 

There are eight other [C"05l/IA07l radio-quiet BL Lac can- 
didates lacking RASS X-ray detections for which we do 
not have deeper Chandra coverage. It is likely some of 
these will also be shown to be low-redshift WLQs upon 
more detaile d X-r ay scrutiny. Many radio-quiet BL Lac 
candidates in lPlOl are probably also WLQs, and we thus 
conclude that the SDSS has discovered a new population 
of low-redshift WLQs. 

We compare our radio-quiet SDSS BL Lac candidates 
to other radio-quiet BL Lac candidates/WLQs in the lit- 
erature. We discuss how high-accretion rates can lead 
to soft ionizing continua and yield spectra with weak or 
absent UV lines but relatively strong Ha and H/3 (but 
those lines may still be weak compared to normal Type 1 
quasars). Near-infrared spectroscopy is required to test 
this properly. Two of the eight objects lacking sensi- 
tive X-ray limits are si milar to the radio-qu iet BL Lac 
candidate discovered bv lLondish et al.l (|2004 2QZJ2154- 
3056). These objects are unlikely to have high mass ac- 
cretion rates (because they do not show H/3), and they 
are too luminous to have weak BELRs because of ex- 
traordinarily low accretion rates. 

The radio-quiet BL Lac candidates should be moni- 
tored for polarization to determine if these objects show 
signs of weakly beamed relativistic jets. Objects with rel- 
atively larger (smaller) a ro (a ox ) values or limits similar 
to those for radio-loud BL Lac objects should be placed 
at high priority for polarimctric monitoring. Similarly, 
searches for flux variability and deeper radio observa- 
tions for objects with poorer a ro constraints would be 
useful. The detection of hot thermal dust emission in 
the infrared would strongly argue against a boosted con- 
tinuum, as beaming would dilute dust emission (e.g., see 
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iDiamond-Stanic et al.| J2009D. 

lLondish et al.l (|2004D successfully model the spectrum 
of 2QZJ 2154—3056 as the combination of an elliptical 
galaxy, a starburst, and a central AGN, and an image 
decomposition shows a host galaxy and AGN. Their op- 
tical image cannot discriminate between an elliptical or 
disk host. Int erestingly, PHL 181 1 appears to live in a 
spiral galaxy () Jenkins et al.|[2005h . BL Lac objects are 
hosted ( perhaps exclusive ly) by giant elliptical galaxies 
(e.g., see lUrrv et al.ll2000h . Thus, host galaxy imaging, 
which could be achieved at least for the lowest-redshift 
radio-quiet BL Lac candidates, might be a useful diag- 
nostic. 

The SDSS has discovered a population of at least 
13 low-redshift (z < 2.2) radio-quiet AGN with weak- 
featured optical spectra, many of wh ich w e believe are 
low-redshift analogs to WLQs. The IP10I optically se- 
lected BL Lac sample adds up to another 71 objects . The 
case of SDSS J0945+1009 (IHrvniewicz et al.ll20loh with 
weak Mg II REW ~ 15 A suggests there are likely even 
more objects yet to be recognized in the SDSS database 
(as SDSS BL Lac searches exclude objects with any emis- 
sion feature larger than 5 A rest-frame) . In all likelihood, 
the SDSS radio-quiet BL Lac candidates are a heteroge- 
neous mix of at least two populations of objects, so it is 
going to take even more observational effort to disentan- 
gle the physical nature of these strange objects. 
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